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Abstract 

Ultracold atomic gases have been used extensively in recent years to realize textbook exam- 
ples of condensed matter phenomena. Recently, phase transitions to ordered structures have 
been predicted for gases of highly excited, "frozen" Rydberg atoms. Such Rydberg crystals are 
a model for dilute metallic solids with tunable lattice parameters, and provide access to a wide 
variety of fundamental phenomena. We investigate theoretically how such structures can be 
created in four distinct cold atomic systems, by using tailored laser-excitation in the presence 
of strong Rydberg-Rydberg interactions. We study in detail the experimental requirements and 
limitations for these systems, and characterize the basic properties of small crystalline Rydberg 
structures in one, two and three dimensions. 

1 Introduction 

Since the development of techniques to cool and trap atomic gases with laser radiation in the 1990's, 
it has been possible to realize and study many-body phenomena in dilute, ultracold gases. Some 
of these phenomena are straight out of solid-state textbooks, such as the Mott-Hubbard transition 
[U[2]: the transition from a "conducting" state to an "insulator" state. In optical lattices, which 
are artificial crystals of light created by standing- wave laser beams, bosonic atoms could be trapped 
with a filling factor of unity, and a superfluid to Mott-insulator transition could be observed. 

One serious shortcoming of such cold atom experiments with regard to the comparison with 
condensed matter physics is that neutral ground-state atoms interact via short-range Van der Waals 
interactions, which are of course much weaker than the Coulomb interactions between electrons in 
solids. A much closer realization of a strongly-coupled condensed matter system would arise if the 
atoms could be separated into ions arranged on a lattice, and surrounding electrons that provide 
the conductance phenomena. 

The question then arises how this could be achieved. Ultracold gases can be made strongly- 
interacting by using Feshbach resonances SJ, however, it does not mean that the system is 
automatically strongly-coupled [5]. These gases are typically very dilute, such that the Van der 
Waals range ro = {mC§/Ti 2 ) l l 4 /2 of the ground-state atoms, where m is the atomic mass and Cq the 
Van der Waals coefficient, is much smaller than the average distance between the particles n -1 / 3 , 
with n the density. The coupling parameter V that indicates the transition between weakly coupled 
and strongly-coupled systems, is given by the ratio of interaction energy E- irA and kinetic energy 
-Ekin- Even for a unitary Fermi gas [6], which is by definition strongly interacting, the coupling 
parameter is smaller than 1: V = \Ei n t/Ep\ ~ 0.6, with E^ n = Ep the Fermi energy. 
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Figure 1: Schematic overview of our experimental set-up, on which the model parameters are 
based, (a) Atoms are cooled and trapped in a Magneto-optical trap, (b) Then an ionization beam 
is utilized to create an ultracold plasma. Due to the accelerator, the electrons and ions can be used 
to create bright beams of electrons or ions. The ultracold particle accelerator ((pictured in (c)) 
provides us with a tool to spatially resolve the Rydberg lattices. 



One way to study strongly interacting and strongly coupled atomic systems is by using ultracold 
Rydberg atoms [7]. Rydberg exitations into electronic s — orbital states interact, similar as the 
ground-state atoms, via Van der Waals interactions. However, the Van der Waals coefficient scales 
very rapidly with the principal quantum number as C§ ~ n , and therefore it is possible to reach the 
regime ro > n -1 / 3 . Robicheaux and Hernandez predicted already in 2005 many-body correlations 
which can rise from the Rydberg blockade effect (also referred to as dipole blockade effect) in a 
disordered gas [HJ, hinting at the possibility of ordered structures. Later, Weimer et al. showed that 
a phase transition could occur from a strongly interacting Rydberg gas to a crystalline phase [9], 
for positive detuning from the ground-state to Rydberg transition. Since then several proposals 
have been made to create correlated systems that involve Rydberg atoms: for instance by means 
of dynamical crystallization of Rydberg atoms starting from ground-state atoms |1Q[ II 1| [T2"]. or by 
weakly dressing ground-state atoms to form super-solid droplet crystals [131 ttH US] ■ 

In this paper, we investigate a very practical scheme which was first proposed by Pohl, Demler 
and Lukin [10], which allows for a laser-assisted, adiabatic self-assembly of a Rydberg atom lattice. 
We investigate the robustness of a tailored excitation scheme to build up a correlated crystal state 
from a disordered atomic gas. Although this paper has a theoretical nature, our investigations are 
based on the present ultracold atom set-up in our group in Eindhoven. The experiment contains also 
an accelerator infrastructure (see Figjl]), which provides us with a tool to detect correlated Rydberg 
crystal structures. The set-up contains a magneto-optical trap for rubidium atoms with appropriate 
780 nm trapping lasers. We routinely trap ~ 10 8 85 Rb atoms at densities up to 10 17 atoms/m 3 
and with temperatures below 1 mK. Both pulsed (10 ns YAG-pumped dye) and continuous-wave, 
single- frequency (300 mW solid-state) 480 nm lasers are present for Rydberg excitation |17j . The 
atom trap is built inside an accelerator, followed by a beam transport system that allows one to 
create magnified images of both ions and electrons onto a microchannel plate and phosphor screen 
assembly that is observed with a cooled CCD camera. The acceleration field can be switched on to 
lOkV/cm on sub nanosecond time scales. This set-up has been used to study electron and ion beams 
generated from the laser-cooled atoms by near-threshold photo-ionization and by field-ionization of 
Rydberg atoms. Both ions and electrons can be detected with a multichannel plate as a measure for 
the total number of Rydberg atoms produced. By directly ionizing the sample, and thus creating 
an ultracold plasma, followed by a rapid extraction of charged particles, we demonstrated that this 
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technique can be used to produce cold bright beams of ions and electrons, with temperatures in 
the beam of the ions of T=lmK and for the electrons of T=10K [17, 18J. 

This paper is organized as follows: In Section [2] we give a theoretical desciption of a cold atomic 
gas in the frozen gas limit, and the corresponding Hamiltonian that describes it. In Section [3] we 
argue how Rydberg crystals can be created by adiabatically chirping the laser frequency. In Section 
[4] we look in detail into four distinct cold atomic systems which we consider for creating Rydberg 
lattices. In Section [5] we present our results for crystal formation using tailored excitation pulses 
for one-dimensional systems. In Section [6] we investigate a few basic crystal properties of two- and 
three-dimensional systems, and finally we conclude in Section [7j 

2 System description 

Throughout this paper we will be working in the frozen gas limit, where the excitation dynamics 
happen on a time scale so fast that we can neglect all motion of the atoms [10|, 119] . i.e., ignore their 
kinetic energy. 

We will investigate four different systems: (i) Atoms in a deep optical lattice, (ii) Condensed 
atoms in a BEC, where the Thomas-Fermi approximation can be made, (iii) Thermal atoms in 
a dipole trap, where temperatures can be in the [iK regime. And finally, (iv) Thermal atoms in 
Magneto- Optical Trap (MOT) conditions, with temperatures typically in the sub-mK regime. 

Under all of the above conditions, the frozen gas approximation is satisfied. Even for the 
highest temperatures in the MOT, it takes tens of /is for atoms to move significantly on the crystal 
scale of 5 — 10/wn |10|. I19j. The Doppler limit, which gives rise to the temperature of the atoms in 
the MOT, can be expressed as a velocity that is about 0.14 m/s for rubidium. With this velocity, 
at the time scales of interest of 1 — 10/xs, we can safely assume that the interaction energy between 
Rydberg atoms is constant. Nevertheless, a dephasing mechanism based on Doppler broadening due 
to the thermal motion of the atoms could potentially frustrate lattice formation. In the following 
sections, we show how crystal creation relies on adiabatic changing of system parameters. In |A"| we 
verify, in a two-particle picture, that Doppler broadening does not affect such dynamics. 

In general, the systems that we describe consist of ./V atoms in a laser field. We assume an 
effective laser coupling between the ground state and a particular Rydberg state, such that we can 
describe each atom as a two level system. The iV-body excitation state can then be represented 
by basis states of the form |n), where n is a iV-dimensional vector with elements m = {0,1}, 
indicating whether the i-th atom is in the ground (0) or excited (1) state. First, we consider the 
special situation in which the N atoms are localised at fixed positions. We will shortly show how 
this picture can be applied to more general systems. The positions of the N atoms are denoted 
with the 3xJV vector r, where rj G R 3 is the position of the i-th atom. In the rotating frame, the 
Hamiltonian governing the dynamics of the excitation state vectors |n) of the atoms is given by 

mm 

N -I N N n 

H exc (r) = - e aA*g + ~ E m OS + 4?) + E ^rAMl (i) 
i=i j=i i=i 1 1 j* 

j>i 

where = |lj) (lj| counts whether the i-th atom is excited, and <j[q = |lj) (0j| and = |0j) (lj| 
change the state of the i-th atom. The first two terms in Eq. ([!]) describe the atom-laser interaction, 
with A the laser detuning and f2 the effective Rabi frequency of the system. The third term equals 
the Van der Waals interaction between atoms in the Rydberg state. 

In practice there are often many individual atoms spaced so closely together, that we can 
safely say that within a certain group of atoms there will never be more than one excitation due 
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to the Rydberg blockade. Moreover, we cannot (e.g. because of detection resolution), or need 
not, distinguish which particular atom carries an excitation. Such groups of atoms can be said 
to coherently share an excitation, and are commonly called superatoms [20]. We can partition 
the space into M volumes \%, k = 1..M, each of which contains a group of M k atoms which 
together represent the £;-th superatom. The Hamiltonian ([!]) can be equivalently applied to describe 
superatoms instead of single atoms, where the positions rj are taken to be weighted averages of 
atomic positions within a superatom, and the Rabi frequency is replaced by a collective Rabi 
frequency n k for the fc-th superatom [S EH E0 1221 E3] : 

n -> n k = y/M^n. (2) 

This collective Rabi frequency reflects the fact that there are multiple atoms in a superatom, and 
that the superatom as a whole is more susceptible to the laser light. This latter fact becomes 
important later on, once we consider dense, inhomogeneous systems. Moreover, in Sec. [5] we shall 
perform simulations of the iV-body system and superatoms will prove to be a useful means for 
reducing the size of the Hilbert space under consideration. 

Having introduced the Hamiltonian for a collection of particles localised at fixed positions, 
it is now time to turn our attention to more general systems such as Bose-Einstein condensates 
with fully coherent wavefunctions, or statistical mixtures found in thermal clouds. To this end, we 
expand the Hilbert space T~L e of excitation configurations |n), with the Hilbert space of iV-body 
spatial wave functions H r , to form the total product Hilbert space Htotai = Hr ® He- 

A convenient basis for this statespace is given by states 

|r) ® |n) , (3) 

where |r) is the state with all N particles localised at positions specified by the 3iV-dimensional 
vector r. The Hamiltonian operator ([!]) is straightforwardly extended to the total Hilbert space: 

H = J d 3JV r|r)£ e *c(r)(r|. (4) 

It is in this final step where the frozen gas approximation becomes apparent: in absence of kinetic 
energy, the Hamiltonian H is diagonal in the spatial components. 

Using the basis states ©, and assuming that all particles start in the ground state |n) = |0), 
we can write the general state of the system at t = as 

|*(0)) = J d 3JV rc (r) |r) <g> |0) . (5) 

When t > 0, this state starts to evolve under the action of the Hamiltonian Q, written down 
formally with the time evolution operator U{t): 

|$(t)) = U{t) |tf(0)) = Jd m v co(r) I7(t)(|r>®|0>). (6) 

The frozen gas Hamiltonian Q however cannot make transitions between states with different 
spatial parts |r) , |r'), with r/r'. Therefore, the time evolution operator of Eq. ^ will only affect 
the excitation part |n) of the many-particle state |r) <g> |n), and not the spatial part |r), so that we 
can write 

I7(t)(|r}®|0>) =\r)®U' xc (t)\0). (7) 
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Here, the time evolution operator U^ xc (t) governs the excitation state of all the particles under 
the action of the Hamiltonian Q, given the fact that they are pinned at definite locations r. The 
general time-dependent state in Eq. Q can then be written as 

|*(t)> = J d 3N rc (r)\r)®\ Xr (t)), (8) 
where we have defined excitation state vectors 

2 N -1 

\Xr(t)) := Ur(t) |0) = Y, C n,r(t) I*) , (9) 

ra=0 

that describe the excitation state of a given atomic configuration |r). The frozen gas Hamiltonian 
Q only affects the excitation part \Xr(t)}, which thus carries the time dependence whereas the 
spatial part |r) remains "frozen". It should be noted that co(r) and c n , r {t) thus satisfy separate 
normalisation conditions given by 

2 N -1 

d? N r |c (r)| 2 = 1, M*)\ 2 = 1- (10) 

n=0 

From the above discussion, we can conclude that the time evolution of a system in the state ([5]) 
can be broken down into the time evolution of the individual, position basis components |r) that 
make up the total state. Each such component evolves completely independent from the others as 
a result of the diagonality of the Hamiltonian Q. Expectation values of observables of interest, O, 
that are independent of momentum, such as the excitation density, are computed for each of the 
individual time evolved configurations. The results are combined in an incoherent sum, in which 
the coefficients co(r) only appear in terms of their modulus squared: 

i3N 



(O) = J d dJV r| Co (r)|^ (r ® Xr (t)| O |r ® X r(t)) . (11) 

Similarly, expectation values of observables in statistical mixtures, such as thermal clouds in a trap, 
are also computed using incoherent sums over individual time-evolved configurations. Here, the 



coefficients |co(r)| in Eq. (11) are replaced by the statistical probabilities of finding a particular 



atomic configuration r. In the remainder of this paper we shall therefore employ Eq. (11) for both 
coherent states as well as thermal clouds, with the understanding that the |co(r)| 2 coefficients are 
statistical probabilities in the latter case. 



3 Adiabatically creating Rydberg crystals 

In the previous Section we have argued how the time evolution of the many-particle state can be 
broken down into the time evolution of the individual, position-basis states |r) ® |Xr(^)) that make 
up the total state of Eq. Q. To see how a crystal state can adiabatically evolve in the frozen 
gas limit, we now single out a particular component configuration |r) of the atoms, and analyse 
the time evolution of the corresponding excitation state |Xr(*))- Since we restrict ourselves to one 
particular configuration |r), we shall omit the spatial ket from the product state notation. 

In the limit of zero laser intensity, O — > 0, the uncoupled states |n) are eigenstates of the 
Hamiltonian ([!]). The energy E n of each state |n) assumes a simple linear dependence on the 
detuning: 

E n = -mA + E int , (12) 
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where m = J2i n i 1S the total number of excited atoms in the state |n), and E{ n t is the total 
interaction energy between the excited atoms in this particular configuration. Clearly, the energies 
of all states with the same number of excited atoms have the same slope as function of the detuning 
A. However, depending on the specific location of the m excitations, each state has a different offset, 
Ei n t, at A = due to the interaction energies. 




A 



Figure 2: (Schematic) Energy of the basis states |r) <g) |n) (thin, gray lines), as a function of the 
detuning and for a given, fixed, choice of the frozen atom positions r. All energies depend linearly on 
the detuning with a slope given by the number of excitations m, and the offset at A = is determined 
by the Van der Waals interaction energy between the Rydberg atoms. For each value ofm there exists 
a specific configuration with the lowest energy (thick, red line), in which the Rydberg excitations 
are typically regularly ordered (schematically pictured configuration on the left). Less ordered states 
have higher energies (example configuration on the right). Note that the atomic positions do not 
change, merely the excitation localisation is different between the two configurations. For non-zero 
laser coupling S7 > 0, the ground state (dashed line) becomes a superposition of basis states and 
separates in energy, which can be adiabatically followd in the direction of the arrow when slowly 
changing the detuning. 

As such, all states with a specific excitation number m form a manifold in the E — A plane, 
as illustrated in Fig. [2] by sets (a "manifold") of gray lines (individual states) with the same slope. 
Now, there are two crucial points to be noted. 

Firstly, the lowest lying state within each manifold is typically an ordered, crystalline state, 
indicated with a solid, red line in Fig. [2} That is, the excitations will tend to be found at regular 
spacings (in Sec. [6] we discuss the crystal structure in more detail). The states with the same 
excitation number, but a less ordered distribution of excitations will typically have a much higher 
energy. This is indicated schematically with two example configurations in Fig. [2] For each value 
of the detuning, there will be a ground state with one specific excitation number, and one specific 
configuration of excitations, which has the lowest energy of all possible states. For — >■ 0, the 
ground state will always be an ordered, crystalline state. 

The second important point concerns what happens when £1 > 0, and regards the possibility 
to adiabatically create a crystal state. As the laser intensity SI is increased from zero, it couples 
the different states |n) through the second term in the Hamiltonian (fiT). Crossings between states 
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from adjacent manifolds in Fig. [2] now turn into anticrossings. In particular, the ground state now 
becomes a superposition of states |n), and separates in energy from the others, as indicated by a 
black, dashed line in Fig. [2j lying below all others. 

A scheme to generate a crystal state would then be as follows. We start the system with a 
large negative detuning, f2 = 0, and all atoms in their atomic ground state, such that we are also 
in the ground state of the many particle system. When the laser coupling is then turned on, and 
subsequently the detuning is raised slowly enough, the many particle state will adiabatically follow 
the ground state |10j . as indicated by the blue arrow in Fig. [2| This process can continue until 
some final detuning is reached, at which point the laser coupling f2 is switched off again. If this last 
step is also performed adiabatically, the system will end up in the O = ground state, which is one 
of the previously discussed ordered crystal states, with a given number of Rydberg excitations each 
of which is localised on exactly one of the frozen atoms, and regularly spaced, forming a crystal. 

The preceding discussion has focused on one particular configuration of atoms |r). However, 
the system state is typically built up from many different configurations, either as a superposition 
of the form <^ in the case of a BEC, or a statistical mixture in the case of a thermal ensemble. 
As explained in the previous section, expectation values of observables for such systems are formed 
by incoherent sums of the results computed for such individual configurations. Typically, crystal 
structures are retained in this summation procedure. 

In the next section we will outline experimental conditions under which adiabatic crystal 
preparation is theoretically possible, followed by Sec. [5] where we will perform simulations of such 
experiments. 

4 Physical systems 

In this section we discuss in more detail the different experimental systems where Rydberg crystal 
creation is possible. We strict ourselves to a (quasi) ID geometry for simplicity as well as for 
computational feasibility, however, it should be noted that the principle of crystal creation is 
straightforwardly extended to higher dimensions. We consider 85 Rb atoms initially in the ground 
state 5s 2 5 , 1 / 2 i which we can excite to the n = 65 Rydberg state ras 2 S , 1 / 2 by using a two-step 
laser excitation scheme. The first step involves an off- resonant 5s-5p 780 nm transition with Rabi 
frequency Q± and detuning hS, and the second step involves a 5p-ns 480 nm transition with Rabi 
frequency 0,2. The two-photon Rabi frequency is then given by Q, = For the n = 65 

Rydberg state, the Van der Waals coefficient is given by Cq = 2.4 10 -58 Jm 6 [21]. We also consider 
a fixed ID excitation volume, of length a=45 fim, which will be defined by the excitation laser 
profiles and / or the spatial dimensions of the atomic system. 

The four main types of systems in ultracold atom experiments are (i) optical lattices, (ii) 
Bose-Einstein condensates, (iii) thermal clouds in an optical dipole trap (high density and low 
temperature), and (iv) thermal clouds in a MOT (low density and relatively high temperature). 
For each of these systems we will show under which conditions Rydberg crystal creation should 
be possible, and we provide some typical numbers as well as investigate possible drawbacks and 
advantages. 

(i) Optical lattices: The conceptually simplest system is formed by ultracold atoms trapped in 
a deep optical lattice. Here we consider a slightly idealised system as an example from which we 
can understand the more complicated systems. We assume that within the excitation volume there 
is a fixed and uniform number of atoms Na per site, and that the excitation volume is sharply 
defined. The lattice is assumed so deep that the atoms are in the Mott insulator phase and there is 
no hopping, such that the atoms are trivially in the frozen gas limit. The system is then represented 
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by a regularly spaced chain of superatoms, each representing Na individual atoms. Typically, Na 
is simply equal to 1, and lattice spacings of 1 /im are commonly realised [25], such that there in 
the order of 40 — 50 atoms within the excitation volume, all experiencing the same coupling to the 
laser field. 

The crystal states in turn are formed by regularly spaced excitations of the atoms. The critical 
detunings A™ +1 for which the lowest energy excitation configuration changes from m to m + 1 
excitations, can be estimated as [TU] 

hAT 1 = § (m 7 - (m - l) 7 ) . (13) 

Here we have assumed that the excitations are regularly spaced, and ignored next-nearest neighbour 
interaction energies. Thus, to create a crystal of m excitations, the final value of the detuning at 
the end of a chirp should lie between A™_ 1 and A™ +1 . 

Finally, we like to mention that crystal creation in optical lattices in a ID geometry has already 
been extensively investigated by Schachenmayer et al. in reference [11], finding interesting effects 
due to descrepancies between lattice and excitation spacing. Similarly, quasi ID systems in a ring 
geometry have been thoroughly studied by Olmos et al. |26| [271 123] who, for instance, make use of 
symmetries to reduce computational costs. 

(ii) Bose Einstein Condensates: The next type of system we consider is a Bose-Einstein 
condensate (BEC) in the Thomas- Fermi limit [2H]. The main difference with the optical lattice is 
a much higher number of particles in the excitation volume, as well as a non-uniform continuous 
density distribution. As we saw in the superatom picture described in Sec.[2j this affects the way the 
system responds to laser light. In the Thomas-Fermi (TF) limit, the density profile n(r) assumes a 

parabolic form as a function of position z along the tube: n{z) = no ^1 — (^) J > where R z is the 

TF-radius of the condensate and no is the central density. The excitation volume can be created 
by a wide 780 nm beam illuminating the entire condensate, and a narrowly focused 480 nm laser 
passing through the cloud, requiring R z = 22.5//m to obtain a 45u m long tube. The diameter of 
this tube will be set to 4 fim, which is reasonable regarding the diffraction limit for this wavelength. 

Typically, no ~ 10 14 cm -3 [201 [29], leading to O(10 4 ) particles in the excitation volume and a 
condensate well within the Thomas- Fermi limit. 

Since the density of the system is not uniform, the many-particle position configurations | r) of 
Eq. Q also have a non- uniform probability distribution cq(t) associated with them. In particular, 
the probability distribution P of the position rz of the leftmost particle is equal to 

P(r L <r)= F(r L ) N , (14) 

where F(r) is the cumulative distribution function of the single particle probability distribution. 



The probability density for the rightmost particle is found by mirroring Eq. (14) around r = 0. 
These probability distributions are important, as the outermost particles of a particular configu- 
ration |r) are certain to carry an excitation in the corresponding ground state. As a result, the 



crystal size a used in calculating the critical detunings A™ +1 in Eq. ( 13 ) will vary between different 
configurations |r), and the critical detunings will likewise assume a probability distribution for their 
values. Fig. (3^a) shows the probability distribution of critical detunings A|, A|, and Ag, relevant 
for creating a 5 or 6 excitation crystal, calculated with 3.6 • 10 4 particles in the excitation volume. 
As opposed to the optical lattice case, the distributions now have a certain width, for instance the 
Ag value has a spread of about 0.5 MHz about its mean value. To create a 6 excitation crystal with 
high fidelity one therefore has to have a final detuning at least 0.5 MHz away from the mean value 
of Ag. The peaks are well enough separated that it is at least always possible to find a suitable 
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Figure 3: Target detuning probability densities for (a) a BEC in the Thomas-Fermi limit, (b) a 
thermal cloud in an optical dipole trap, and (c) a thermal cloud in a magneto- optical trap. The 
critical values for the optical lattice, as given by Eq. (13) are superimposed as vertical, dashed lines. 
Insets: probability distribution for the locations of excitations in the 5-excitation ground state. 



detuning at which we are certain to have a definite excitation number in the ground state. As 
we will shortly see in the case of thermal clouds, this need not always be the case. The inset of 
Fig. [3^a) shows the probability distribution for the locations of the excitations in the 5-excitation 
ground state. Due to the probability distribution of the outermost particles, the peaks are not 
infinitely sharp but have a finite width. 

(Hi - iv) Thermal clouds in optical dipole traps and MOTs: Thermal clouds with Gaussian 
density distributions form the starting point for many experiments in ultracold atoms. There is a 
great variety in sizes and density, so we will treat two specific cases, clouds confined in an optical 
dipole trap (ODT) and clouds contained in a magneto-optical traps (MOT). 

Clouds confined in an ODT are much smaller in size than those in a MOT, and they reach 
much higher densities. Typical peak densities no in ODTs range between 10 12 - 10 15 cm 3 [201 130] . 

We consider a cylindrically symmetric density of the form 

n(r ' Z) = (27r)I 2 a3 ex P(" r2 /2^V 2 2 - z 2 /2a 2 z ), (15) 

where k is the aspect ratio of the cloud, N the total number of atoms in the cloud, and a z the width 
of the distribution in the z-direction. To create an excitation volume of 45//m long in the ODT, 
we can proceed in the same fashion as was the case of a BEC, by using a defocused 780 nm beam 
illuminating the entire sample, and a narrowly focused 480 nm beam aimed through the center of 
the cloud along the z-axis. For typical numbers iV = 4 • 10 5 atoms, an aspect ratio of k = 10 and a 
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width of a z = 7/xm, we get 2000 particles in a tube of 4//m in diameter. These values are precisely 
such that the expectation values of the positions of the outermost particles in the tube are spaced 
45/um apart. 

As was the case for the BEC, the nonuniform probability distribution for the particle density 
causes the critical detunings of Eq. (13) to assume a probability distribution as well. Fig. |3^b) 
shows the probability density for the critical detunings A|, A|, and Ag relevant for creating a 5 or 
6 excitation crystal in our excitation volume in the optical dipole trap. Apparently, in this case 
the probability densities overlap, which means that there exist no values of A in which the ground 
state is certain to have one particular number of excitations, there is always a finite probability to 
have one more, or one less excitation than expected. The inset of Fig. [3^b) shows the probability 
distribution for the locations of the excitations in the 5-excitation ground state. 

Clouds in a magneto-optical trap (MOT) are typically much larger in size, with diameters of 
0(1 mm). The 45/im excitation volume therefore needs to be created by shaping the intensity 
profiles of the excitation lasers. By passing the 780 nm laser through an aperture, it is possible to 
create a beam of width 45/im and approximately uniform intensity locally within the atom cloud. 
Intersecting this beam with a small, focused beam with diameter 4//m of the 480 nm laser near the 
center of the atom cloud results in the desired ID excitation volume. Inside this volume the atomic 
density is practically uniform, and with typical densities of a MOT of n = O(10 10 — 10 11 )cm -3 there 
are approximately 10 to 60 atoms at uniformly distributed positions contained within the excitation 
volume. Figure [3]^c) shows the critical detunings for such a system for the particular case of 30 
atoms. The probability distributions form peaks with distinctly asymmetric shapes, which are a 
result of the sharp boundaries of the excitation volume. The inset of Fig. [3^c) shows the probability 
distribution for the locations of the excitations in the 5-excitation ground state, which also exhibit 
the characteristic asymmetry due to the sharp boundary of the excitation volume. 

As a final note to this section, we would like to mention that it is in principle possible to 
change system parameters such as the dimensions of the excitation volume. However, due to the 
6-th power scaling of the Van der Waals interaction energies this has a large effect on the other 
system parameters. For instance, changing all length scales by a factor A is equivalent to changing 
all time scales by a factor of A 6 , and consequently all frequencies such as f2 and A need to be scaled 
by a factor A~ 6 . Even a simple change of the excitation volume length from 45/mi to 90/im would 
therefore require all timescales to be multiplied by a factor of 64 in order to obtain the same results. 
The chirps that we will construct for our particular excitation volume and particular experimental 
setup can therefore be said to be truly tailor made. 



5 Performing the chirp 
5.1 Chirp properties 

Having outlined the properties of the physical systems of interest, we proceed by composing a 
realistic chirp of the laser detuning A and intensity 0. The basic ingredients for such a chirp are 
as follows. 



1. Start by switching on the laser and increase the coupling strength £1 from zero up to some 
final value. This is done at a large, negative detuning with all particles starting in their 
atomic ground state, coinciding with the many body ground state. 

2. Once the coupling £1 is sufficiently high, the detuning can be increased from its initial value 
towards its final target value Af between A™_ x and A™ +1 . 
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3. Switch off the laser, by lowering the coupling strength back to £1 = again. 

The above three steps should be carried out as adiabatically as possible, yielding a crystal state 
at the end of the chirp, such as depicted in the insets of Fig. [3j In Figure |4j panels (a), we 
show various typical chirps, modified by a computer algorithm to maximise adiabaticity given a 
constrained time frame of ~ 5//s in which it is to be carried out. Each chirp is tailor made for a 
specific system. The remaining panels (b) and (c) of Fig. [4] will be discussed in Sec. |5.3[ 

Starting with step (i) of the chirp, the first feature to be noted is the shape of Q(t) in the 
first section, which is roughly exponential. As the coupling is still small, the energy levels are 
closely spaced together, and only small changes to the Hamiltonian can be made without making 
transitions to excited states. When the coupling is stronger, the ground state gets pushed to lower 
energies and separates further from the other energies, and as a consequence we can make larger 
changes to the Hamiltonian and increase the rate of change of the coupling. As the final value is 
reached, the ramp abruptly stops and Q remains constant. Since the Schrodinger equation is a first 
order differential equation in time, such sharp changes in slope do not affect adiabaticity and there 
is no need to smoothen the chirp. 

When in the second step the detuning is increased, we found that the maximum slope varies 
non-monotonically throughout the chirp. Regions where the detuning can only be changed relatively 
slowly are not necessarily correlated to the critical detunings A™ +1 described in Sec. § This is 
a sign of a complicated underlying energy landscape, whose details are beyond the scope of this 
paper. 

The most critical step of the chirp is the final step, switching off the laser and lowering Q back to 
zero again. When the detuning is large and positive, the energy scales in the system become larger, 
such that the coupling strength becomes relatively smaller. Moreover, in the case of a non-uniform 
density, such as the ODT or BEC, the y/Mk enhancement of the coupling in Eq. Q necessitates 
a very fine control over the final stages of switching off the coupling, as even small intensities can 
still couple strongly to the central, high density region of the excitation volume. Because of these 
two phenomena, step (Hi) takes up by far the greatest amount of time of the total crystal creation 
procedure. Typically, the time scale for switching off the laser is an order of magnitude larger than 
the combined time for switching it on and chirping the detuning. Again it should be noted that 
the rate of change in 0, varies, as another sign of a complicated underlying energy landscape. We 
also found that it is favourable to decrease the detuning slightly while ramping down the intensity. 

5.2 Simulating excitation dynamics 

We now have all the necessary ingredients to perform a realistic chirp and simulate an attempt to 
create a Rydberg crystal in each of the four systems described in Sec. [4j an optical lattice, a BEC 
and thermal clouds in an ODT and a MOT. We aim to create 5 or 6 excitations in total, due to 
limitations imposed by the computational complexity. 

Following the ideas outlined in Sec. [2j the simulation can be broken down into two parts. The 
first part consists of simulating the time evolution of the excitation dynamics of a given set of atoms 
at fixed positions. Then, in the second part a Monte Carlo integration over particle distributions 
is performed, summing many simulations for different atomic positions weighted by the probability 
that a certain configuration occurs. 

For the simulation of the excitation dynamics for a fixed configuration, we perform a simulation 
of the full quantum many particle state, governed by the Schrodinger equation with Hamiltonian 
([!]). We use a fifth order adaptive stepsize Runge Kutta time integration [31J, with the Cash- 
Karp parameter scheme [32J. For the optical lattice we assume 45 particles in the excitation 
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volume, whereas in the case of the BEC and the ODT we assume N = 0(1O 3 - 10 4 ) atoms in the 
excitation volume, and 30 in that of the MOT. Since the parameter space for N two level systems 
scales as 2^, we need to make some approximations to make the simulations tractable. The first 
approximation consists of reducing the O(10 3 — 10 4 ) atoms in the ODT and BEC to 40 superatoms, 
by employing the algorithm suggested in Ref. [8] where we recursively replace the closest spaced 
pair of (super)atoms by a weighted average position. In the end, we replace the couplings to the 
superatoms by the ^/M^VL coupling of Eq. Q , where denotes the number of atoms represented 
by the &:-th superatom. But even the resulting 2 40 states space is too large to handle, and we 
reduce the number of states further following the suggestions put forth in Refs. [9l [33l I10j . where 
we omit all states with more than a certain maximum number of excitations (> 8 in our case), and 
also omit all states where the excitations are spaced too close together. This latter approximation 
is made under the assumption that these states will never be significantly populated due to the 
Rydberg blockade. These three approximations reduced the state space to O(10 5 ) states, a number 
low enough to allow for a few hundred simulations per day on an 8 core system. 



5.3 Results 

For each of the systems described in Sec. [4] we have performed simulations of a tailor made chirp 
of laser detuning and intensity. The laser chirp is pre-determined per system by an algorithm that 
attempts to maximize adiabaticity, while keeping constraints on the total duration of the chirp. In 
this way, we have constructed chirps that create 6 excitations in the optical lattice and BEC, and 
5 excitations in the thermal clouds in an ODT or MOT. All chirps have a duration of about 5 /is. 

For the optical lattice, the simulation only needed to be performed for a single configuration 
of atoms, since we assumed that the ground state atoms are at very well localised positions. For 
the other three systems we have performed many Monte Carlo simulations to account for the 
probabilistic distribution of the ground state particles. At the beginning of each simulation, we 
draw new particle positions from the appropriate particle distribution, and then simulate the time 
dynamics when this particular atom configuration is subjected to the laser chirp. In total, we 
have performed 600 Monte Carlo samplings for the BEC, 1650 samplings for the ODT and 3300 
samplings for the MOT. 

After all the simulations are completed, the results can be aggregated in order to compute 
expectation values of various observables. In particular, for each system we compute expectation 
values for the number of excitations, the excitation density, and the pair correlation function. 

These expectation values can be constructed using the excitation number operator ijj\{r)ijj e (r), 
which records the number of excitations measured at the (ID) distance r from the origin (not to 
be confused with the many particle position vector r). We define this operator through its matrix 
elements in the basis ([3]): 

(r ® n'| $(r)$ e (r) |r ® n) = 5{r - r')<Wl r (r, n), (16) 

where l r (r, n) is an indicator function, which takes on the value of 1 if the state |r) £g> |n) has an 
excited particle at position r, and otherwise. The corresponding expectation value is found from 



Eq. 

2 N -1 



<$(r)&(r)> = f d m r £ |c (r)| 2 \c n , r (t)\ 2 l r (r,n). (17) 

J 71=0 

Note that the coefficients co(r) follow from the Monte Carlo sampling, whereas the coefficients 
Cn,r(t) are calculated in the simulation of the time dynamics. 
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Using the operator ipl(r)tp e (r) we can additionally compute the second-order density correlation 
function [34J 

G^(r, r') = (^i(r)Mr)^t(r')Mr')} - 8(r - r'){^ e (r)^ e {r)) , (18) 

which measures the conditional probability of finding an excitation at point r' , given that there 
already is an excitation found at r. It is more convenient to compute the integrated second-order 
density correlation function 

G^{r) = J G^(r,r + R) dR, (19) 

which measures just the probability of finding an excitation a distance r away from any other 
excitation. Calculated on a discrete grid, the correlation function can be thought of as a 
histogram of pair separations. 

The final observable we calculate is the probability p(m) of finding m excitations, which can 
be calculated as 

2-^-1 

p(m) = / d 3iY r|c (r)| 2 £ |c n , r (i)| 2 l m (n), (20) 

where we have defined an additional indicator function l m (n), which takes on the value 1 if the 
state n has m excitations, and otherwise. 

We are now ready to discuss the results of the simulations. Figure [4] shows the calculated 
expectation values for each of the four systems we are interested in: (i) optical lattice, (ii) Bose- 
Einstein condensate, (iii) thermal cloud in an ODT, (iv) thermal cloud in a MOT. For each system, 
the top panel, (a), shows the chirp of the laser intensity f2 and detuning A, the second panel, labeled 



(b), shows the expectation values p(m) of Eq. (20) of the number of excitations, the third panel 



(c), shows the excitation density P7t , where all expectation values are plotted as a function of 
time. 

There are a number of similarities and differences to be noted between the four systems. 
Starting with the top panel (a), we see that the shape of each chirp follows roughly the three- 
step scheme described in Sec. |5.1[ although the length of each segment varies between systems. 
Typically, the third step of switching the laser off takes the longest time. 

Panels (b) and (c) of Fig. [4] show clearly what happens during the three steps of the chirp. 
When the laser is switched on in the first step, the initial state |0) with all atoms in the ground 
state becomes a superposition of states with or 1 excitations. Then, as the detuning is increased 
while the intensity is at its maximum, states with higher numbers of excitations start to mix into 
the system state. No particular number of excitations has a high probability of occurring (panel 
(b)), and the location of the excitations is also spread out across the excitation volume (panel (c)). 
Only when in the final step the laser is switched off do the excitations localise at regularly spaced 
positions, and the Rydberg crystal takes shape. First, the outermost excitations localise, and then 
excitations crystallise progressively inwards towards the center of the excitation volume. In panel 
(i.c) for the excitation density of the optical lattice, one can distinguish the individual atoms at 
each lattice site, and see how the excitations are initially shared among all atoms and crystallise 
only at the end of the chirp. 

The middle section of the chirp takes the shortest time in the systems with high particle 
numbers, the BEC (ii) and the ODT (iii), as a result of the enhanced y/Mk coupling in the superatom 
picture [see Eq. Q]. The effect of this enhanced coupling is also clearly visible in panels (a) of Fig. 
[4] for systems (ii) and (iii), where a disproportionally long time is needed to ramp down the last 
few percents of laser intensity. The center panels (b) show that in this long flat tail of the intensity 
chirp there are still large changes in the excitation number probabilities. 
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Figure 4: Simulation results for Rydberg crystal creation in each of the four systems discussed in 
the text: (i) optical lattice, (ii) Bose-Einstein condensate, (hi) thermal cloud in an ODT, and (iv) 
thermal cloud in a MOT. All results are shown as a function of time. In each subhgure, the top 
panel, labeled (a) shows the laser intensity O and detuning A, the center panel, labeled (b) shows 
the probability p(m) of finding m excitations. The bottom panel, (c), shows the spatial excitation 
density (^J(r)^ e (r)) during the chirp, indicated with colors ranging from white (lightest) to blue 
to red (darkest), where darker colors indicate a higher probability of finding an excitation at that 
position and time. 
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Figure 5: (a): Excitation density, and (b): correlation function G^ 2 \r) at the end of the chirp for 
each of the four systems under consideration: (i) optical lattice, (ii) Bose-Einstein condensate, (iii) 
thermal cloud in an ODT and (iv) thermal cloud in a MOT. Note that the optical lattice and BEC 
have 6 excitations, whereas the ODT and MOT have 5 excitations. 



Fig. |5] s hows the excitation density (17) at the end of the chirp (column (a)), and correlation 
function (19) at the end of the chirp (column (b)), for each of the four systems (i) - (iv) of interest. 

Starting with the excitation density in the optical lattice in the top left panel (a.i), we see that 
the excitations are more or less perfectly localised on individual lattice sites. Upon close inspection 
there are small probabilities discernable, for finding the excitation at neighbouring lattice sites, 
indicating that the time evolution has not been perfectly adiabatic. The excitation density of the 
BEC in panel (a.ii) also exhibits distinct and well-separated peaks. However, the excitation density 
of the thermal cloud in an ODT in panel (a. iii) shows peaks that, although clearly visible, are not 
entirely separated spatially. This is due to the boundaries of the excitation volume not being as 
sharply defined as those of the BEC. In the ODT, the spread in the position of the outermost 
particles is largest of all four systems. Even the MOT excitation density in panel (a.iv) shows 
better separated peaks. Due to the sharp boundary of the excitation volume in this system (see 
Sec. [4| the peaks show a roughly similar shape to those of the BEC. However, they exhibit a larger 
spread around their mean positions, reflecting the much lower number (30) of particles present in 
the excitation volume. In effect, the low particle number offers fewer positions for the excitations 
to localise on, and the optimal positions leading to the lowest energies are often not occupied in a 
particular Monte Carlo simulation. 

Finally, moving to column (b) of Fig. [5j showing the correlation function for each system, we 
can extract the following information. Interpreting the correlation functions G^ as histograms 
for the distances between pairs of excitations, we clearly see that there are never two excitations 
within one peak of the excitation densities in column (a). Should this have been the case, then we 
would find nonzero probabilities at distances smaller than the inter-peak distance of ~ lOfim. 

The first peak in appearing around 10/um is the distribution of pair separations of exci- 
tations in neighbouring peaks in the excitation density. Likewise the second peak in G^ is the 
distribution of pairs in peaks separated by a single other peak, followed by the pair distributions of 
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excitations in peaks separated by two peaks, etc. The height of the peaks in the correlation function 
directly reflects the number of pairs of peaks separated by 0, 1, 2, .. other peaks, leading to a 
linear decay in height. The width of the peaks signifies how regularly the excitations are spaced. 
Again, in case (iii) of the ODT the peaks are the least sharply defined, as a direct result of the 
excitation volume not being bounded sharply. In case (i), the optical lattice, there are small sub- 
peaks visible again, which are correlations between the large peaks and sub-peaks of the excitation 
density in panel (a.i). 



6 Higher dimensional systems 




Figure 6: Examples of 2D crystal structures as function of Rydberg atom number m. Initially, 
when there are few Rydberg atoms, all excitations reside on the boundary of the excitation volume 
in order to maximise the interparticle distance. As the number of Rydberg atoms confined in the 
circular volume increases, the inner region of the excitation volume starts to fill as well. At m = 7 
the first particle appears in the center, followed at m = 10 and m = 12 by a second and third 
one respectively. At m = 18, an additional shell starts to form at the center, progressively filling 
up (m = 24,26,..) until there are so many particles in the interior that boundary effects become 
negligible and a transition to bulk behaviour is observed, where a hexagonal lattice is formed in 
the central region (m = 61 and m = 129). The corresponding energies of these configurations are 
shown in Fig. (9^a). 



In the previous section, we focused in detail on the adiabatic creation of one-dimensional 
Rydberg crystals. However, the dimensionality of the crystal structures can be controlled by the 
shape of the volume set up by the excitation lasers. As ID structures are produced by a needle- 
shaped light field, so can 2D structures be created by a sheet of light intersecting the atomic cloud. 
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3D structures can arise in a light volume where excitations are not prohibited by dimensional 
restrictions. In this section we take a closer look at what type of 2D and 3D structures may 
arise when applying the same type of laser excitation techniques in 2D and 3D light volumes, and 
investigate the energy of the ground states as function of the number of Rydberg atoms. This 
latter quantity is of importance when determining the critical detunings A™ +1 described in Sec H 
at which the system transitions from an m-excitation ground state to an m + 1-excitation ground 
state. 



Figure 7: Example of a 3D crystal structure, (a) displays the structure as a whole, while in (b) the 
inner shells can be seen, with an indication of the average distance of the shell to the center of the 
structure. 

While in the ID case the crystal structures are basically similar to regularly spaced chains 
of Rydberg excitations, in higher dimensions the structures can be quite complicated. However, 
for simplicity, we assume that the crystal structures with m excitations are confined by a circular 
plane (2D) or spherical volume (3D), which have hard wall boundaries and where the interior can 
be considered as homogeneous with no external potential. 

We calculate the optimal crystal structures by using a Conjugate Gradients |31j minimization of 
the energy of m point particles interacting via a repulsive two-body Van der Waals Ce/r 6 potential: 



The minimization is performed with respect to the particle positions r,. Since this energy represents 
a very complex 2m (2D) or 3m (3D) dimensional energy landscape, with many local minima, we 
have performed hundreds of minimizations for each Rydberg particle number m with random initial 
conditions, and recorded the lowest energy configuration thus found as a good approximation for 
the global minimum. There exists a large body of literature in the field of Coulomb crystals where 
similar procedures are found to find crystal structures of trapped ions, see e.g. Refs. [351 133 113 
EHIEH], and references therein. 

In Figs. [6] and [7] examples can be seen of crystal structures of Rydberg excitations, for a given 
number of m. Typically, we observe in 2D that at low particle numbers m all particles are located 
on the boundary of the volume, maximizing the distances between the nearest neighbours. As 
m is increased the boundary fills up, until at a certain point it becomes energetically favourable 




r 
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Figure 8: Shell formation for a 2D (left) and 3D (right) Rydberg lattice. The Rydberg atom number 
density is indicated as a function of the radial distance to the center of the system (horizontal axis) 
and for a varying total number of Rydberg atoms (vertical axis). Each horizontal slice of the 
figure corresponds to particular configuration, as indicated for instance in Figs. [6]and[7j When the 
number of atoms is increased, additional shells appear (blue / red fuzzy lines) that move toward 
greater radial distance to make room for yet more inner shells. A large part of the atoms resides 
on the boundary of the system, corresponding to the bright line at radial distance r = 1. In the 
2D system, the onset of lattice formation and breakdown of shell structure can be observed near 
r = for high atom numbers where the shell structure becomes more diffuse. 



for the new particle to be placed in the center of the volume. As the particle number is increased 
further, more particles will be located near the center, forming a ring shaped structure. This process 
continues until, in turn, a particle appears in the center of this ring again, starting the cycle anew. 
For even larger particle numbers, there are so many particles located in the interior that boundary 
effects become negligible and the particles assume a hexagonal lattice structure typical of bulk 
behaviour. 

Generally speaking, for larger numbers of excited atoms a competition occurs between lattice 
structures in the interior, and shell structures near the boundaries of the system, an effect familiar 
from Coulomb crystal studies |35t [38] ■ Positions on the boundary of the system are energetically 
favourable due to the (partial) absence of neighboring atoms, leading to boundary effects (shell 
structures), propagating inward until the boundary is sufficiently far away such that bulk behavior 
(lattice structures) start to dominate. For a sufficiently large number of Rydberg atoms, the 
structure can then be regarded as a (quasi) periodic, homogeneous lattice. However, for small 
atom numbers we enter the regime of mesoscopic physics, where we create symmetric few-body 
clusters. 

In 3D systems we see the same chain of events, with consecutive shells forming when the 
particle number is increased. Here, we expect lattice formation in the interior only at much higher 
Rydberg atom numbers since surfaces of spheres can accomodate much more particles than the 
circular boundary in the 2D case. We have verified that up to m = 1000 we still see distinct shell 
structures in 3-D. 

As a function of the number of Rydberg atoms in the lattice, we show in Fig. [8] where the 
atoms will be located as a function of the radial distance to the center of the excitation volume. 
Both in 2D and 3D the formation of shells is visible for increasing Rydberg atom number m. 
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Figure 9: Minimum energy E as function of the number of Rydberg atoms m, for the distribution 
of m atoms in (a) a two-dimensional circular plane, (b) a three-dimensional spherical volume. The 
insets display the behavior for smaller atom numbers. The continuous lines represent fits according 
to the mesoscopic (blue, solid line) and bulk (red, dashed line) scaling laws. 



In Fig. [9] the minimum energy E(m) is shown as a function of the particle number m, for the 
2D and 3D case. The energy solutions show three distinct regimes: the few-particle regime where 
boundary effects dominate, the bulk regime with many particles, and inbetween the mesoscopic 
regime which connects the two limits. Each regime has a different scaling of energy versus particle 
number. In the bulk regime the energy scales as m 4 in 2D and as m 3 in 3D. This scaling behavior 
assumes a homogeneous distribution of the particles with a mean interparticle separation, and where 
each particle has a fixed number of neighbours. Furthermore, we assume only nearest neighbour 
interactions. In the few-particle regime, edge effects dominate, which results in a relatively larger 
interparticle separation as in the bulk regime. However, despite a lower than bulk energy, the 
energy increases more rapidly with m. In the mesoscopic regime, the system is inbetween the edge- 
dominated and bulk dominated regimes, which results in energy scaling laws as ^/mm 4 for 2D and 
-y/mm 3 for 3D. 

7 Conclusions and outlook 

While crystalline phases have been predicted for strongly interacting Rydberg gases, it is not 
straightforward to achieve such an ordered structure starting from a disordered gas. A practical 
scheme to enter this regime has been proposed by Pohl et al. [TO], using chirped laser pulses. In 
this paper, we investigate the feasibility to create one-dimensional Rydberg crystals in four distinct 
cold atomic systems, which vary in density and temperature, and demonstrate the importance 
of an excitation scheme which is exactly tailored to the specific system parameters. Even for a 
system that is relatively low in density and high in temperature, such as a magneto-optical trap, 
it is possible to design an excitation scheme that creates spatially ordered structures. A successful 
experimental realization of one-dimensional Rydberg crystals allows for a direct comparison to 
predicted analytical properties in a lattice [10]. For crystalline Rydberg structures in two or three 
dimensions, we calculate some basic properties such as the spatial structure, the corresponding 
ground state energy, and shell structure formation as function of the Rydberg atom number. 
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A A closer look at the frozen gas approximation 

Throughout this paper, we assumed that the system is sufficiently cold such that we can make the 
frozen gas approximation, and hence neglect the kinetic energies of the atoms during the course 
of a typical experiment lasting a few microseconds. In a typical experiment, the laser detuning is 
swept through a few MHz within a few fj,s, and in this Appendix we investigate whether thermal 
motion under MOT conditions can significantly affect the final state after a sweep. 





Figure 10: Velocity-averaged two-particle squared amplitudes as function of time, during the course 
of an adiabatic chirped transition from the ground state to the excited Rydberg state. Shown are 
\c gg \ 2 (red), | c 9e 1 2 + c eff 1 2 (blue), and c ee | 2 (gray). From t = — 4/is to t = 4 /us, 5 is chirped from -15 
MHz to 15 MHz in a linear fashion. Q is held constant at 1.4 MHz. The interaction energy Vi n t 
of the c ee state is (a) 0.2 MHz and (b) 7.8 MHz, corresponding to a distance of about 11 and 6 
/im respectively. The absolute squares of the amplitudes are averaged over a Maxwell-Boltzmann 
distribution with a temperature corresponding to a Doppler limit of 0.5 MHz. Also shown are the 
instantaneous eigenstates of the ideal system with v\ = v 2 = (black line). The decaying Rabi 
oscillations can be seen, which are found to damp out in a Gaussian fashion with a time constant 

Note that the oscillations in (a) experience an adiabatic 



r ~ 0.45 //s, in agreement with Eq. 26 
revival at t 



4/is because of the symmetry of the sweep. 



We consider a two particle two-level system with a time-dependent Rabi frequency and detun- 
ing, described by the coupled set of optical Bloch equations for the coupled two-particle amplitudes 
l c ffff| 2 ' l c se| 2 , |c e9 | 2 , and |c ee | 2 , where subscripts g and e denote ground or excited states. The 
time-dependent equations are given by 



.dc§^) = (/)+<,,(/)). (22) 



ACgejt) Q(t) 

5 dt "2 



(c gg (t) + c ee (t)) + (-A(t) + kv 2 )c ge (t), (23) 
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.dc^t) = 9^^ Cgg (t) + Cee (t)) + {-A(t) + kv l )c^(t), (2 1) 



i dC ^ = ^(c 9e (t)+c efl (i)) + (-2A(t) + ^ 1 + ^2 + ^)c ee (t). (25) 

The Doppler shifts kvi are added to the detuning, which allow for different velocities for particle 
i = 1,2. A conservative estimate for these shifts is of order 3 • 10 6 rad/sec or 0.5 MHz, resulting 
from the effective transition wavelength of 300 nm, that gives rise to the laser wavenumber k. The 
dephasing that arises is a result of the distribution of velocities vi for an ensemble of atoms, which 
we assume to obey Maxwell-Boltzmann statistics, and this will damp Rabi oscillations. The time 
scale associated with this damping is 

T = V tfk*k B T ' (26) 

with m the atomic mass and ks the Boltzmann constant. However, damping will not be a limiting 
factor for the experimental conditions which we consider in this paper, provided that the chirp is 



wide enough, as can be seen in Fig. 10 . Here an adiabatic transition is shown for two ground state 
atoms to their excited Rydberg states, by presenting the velocity averaged two-particle amplitudes. 

c 9a\ 



The calculation is done for two different interparticle separations, and the squared amplitudes |c„ '*-' 



\c ge \ 2 etc. are averaged over Maxwell-Boltzmann distributions, corresponding to the Doppler tem- 
perature, for velocities V{. The figure shows damped Rabi oscillations. Nevertheless, the amplitudes 
are able to follow adiabatically the instantaneous eigenstates if the chirp is slow enough, and wide 
enough, meaning that |A(ii n i t i a i) — kv\ 3> f2 (initial) and |A(ig na i) — kv\ 3> fi(ifmai)- in conclusion, 
since the Rydberg crystal creation relies on adiabatic processes, dephasing of Rabi oscillations does 
not inhibit the crystallization process. 
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